www.gusucode.com > 循环自相关函数工具箱源码程序 > matlab代做 修改 程序循环自相关函数工具箱/cyclostationary_toolbox/ex3.m
clear; clc; t=0:1/600:1023/600; A=1;B=1;fb=7;fa=100; x=A*(1+B*cos(2*pi*fb*t)).*cos(2*pi*fa*t); figure(1) plot(t,x); axis([0 1 -2 2]); %alpha=2*pi/300; T=300; max_tau=0; R=cyclic_autocorrelation_fast(x,T,max_tau); y=R.*conj(R); % ff=600*(0:1023)/512; ff=0:1:300-1; figure(2) % plot(ff(1:150)/(2*pi),y(1:150)); plot(ff,y);